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We consider the behaviour of a quintessence field during an inflationary epoch, in order to learn 
how inflation influences the likely initial conditions for quintessence. We use the stochastic inflation 
formalism to study quantum fluctuations induced in the quintessence field during the early stages of 
inflation, and conclude that these drive its mean to large values (> O.lmpi). Consequently we find 
that tracker behaviour typically starts at low redshift, long after nucleosynthesis and most likely 
also after decoupling. 
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I. INTRODUCTION 



In order to fit the impressive array of observational 
data now available, it is widely believed that a model of 
the Universe should feature a period of inflation in its 
early history, leading to the generation of density pertur- 
bations via quantum fluctuations, and that the present 
expansion rate should be accelerating. While the latter 
phenomenon is often considered due to a cosmological 
constant, it is at least as attractive to presume that the 
acceleration is driven by the same mechanism usually ex- 
ploited to give early Universe inflation, namely domina- 
tion by the potential energy of a scalar field. Such models 
are known as quintessence models. 

An important class of quintessence models are known 
as tracking models [1, 2, 3, 4, 5], where the late-time evo- 
lution of the field has an attractor behaviour rendering it 
fairly independent of initial conditions. However, despite 
the existence of tracking behaviour the details of the ini- 
tial conditions may yet be important. For example, in 
several models of quintessence the tracking solution can- 
not be achieved too early; big bang nucleosynthesis is 
spoiled if the quintessence field has too large a density 
at that time [3, 6]. More obviously, a quintessence sce- 
nario will not work unless the initial energy density is at 
least as large as is required by the present. It is there- 
fore useful to have further guidance as to the likely initial 
conditions for quintessence. 

The simplest assumption concerning quintessence is 
that it is a fundamental scalar field (rather than a low- 
energy composite field), and as such was already present 
during the inflationary epoch in the early Universe. If 
the quintessence field is sufficiently weakly coupled that 
it is not affected by the inflaton decays ending inflation, 
its possible initial conditions are restricted by its dynam- 
ics during the inflationary period. Note that we are not 
considering the situation where the quintessence field and 
the inflaton are the same field [7]. 

In this paper we carry out a comprehensive study of 
the influence of inflationary dynamics on the quintessence 
field. We choose parameters so that the quintessence sec- 
tor matches the observed acceleration of the Universe, 
and the inflaton sector generates suitable perturbations 
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to initiate structure formation. We investigate both clas- 
sical and quantum dynamics of the quintessence. We 
assume a flat Universe throughout. 



II. MODELS AND OBSERVATIONAL 
CONSTRAINTS 

A. Quintessence 

The quintessence model is defined by a scalar field Q 
evolving in a potential V(Q). Constraints on such mod- 
els from observation have recently been considered by 
various authors [8, 9]; the precise details are not impor- 
tant for our considerations but for definiteness we use 
the results of Bean et al. [9], who recently combined con- 
straints on such models from type la supernovae, CMB 
peak positions and large-scale structure surveys. They 
found one-sigma constraints on the quintessence density 
parameter and its pressure-density ratio of 

0.57 < Q < 0.72 and - 1 < w Q < -0.85 . (1) 

With this, and assuming the Hubble constant to be given 
by h = 0.72 ± 0.08 [10], it is possible to find constraints 
on the parameters of the quintessence models. To lo- 
cate viable regions of parameter space, we approximate 
these three constraints as Gaussian distributed (in the 
case of wq a half-Gaussian centred at wq = — 1) and in- 
dependent. This defines a confidence region in the space 
of parameters, and using numerical simulations we can 
translate the confidence region into other variables such 
as the density of non-relativistic matter today and the 
parameters of the quintessence potential. 

Our main discussion will focus on quintessence mod- 
els with inverse power-law potentials, as originally intro- 
duced by Ratra and Peebles [2] 



V(Q) = V 



\ mpi 



(2) 



In order to simplify some later expressions we define the 
notation X = X/mpi, where X is any field or variable. 
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FIG. 1: The 68% and 95% confidence regions for an inverse 
power-law model of quintessence. 



Using the procedure described above and projecting the 
result onto the (3—Vo plane gives the allowed region shown 
in Fig. 1. We see that /? < 1 and V ~ 5 x lCT 124 ™^ at 
68% confidence, with a broader region at 95% confidence. 
In defining our models, we allow (3 to vary and for each 
choice fix Vq at its best-fit value. 



B. Inflation 

We take inflation to be driven by a scalar field (f> (the 
inflaton) in a potential W((j)). Because of the small value 
of Vq , the inflaton dominates the evolution and therefore 
we can use the standard constraints on its parameters. 
The inflation model is to be chosen so as to reproduce 
the amplitude of perturbations needed to generate the 
observed structures, which requires [11] 
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(3) 



This is to be evaluated at ^ bs, the value of the field cor- 
responding to the epoch when the scale of the observable 
Universe crossed the horizon. 

In this paper we will consider potentials where inflation 
ends by violation of the slow-roll conditions, so the slow- 
roll approximation 



(&W/d<j>X 
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< 1 



n{<t>) 
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gives the end of inflation and the number of e-foldings 



(4) 
(5) 



N ~ 



8tt 



pi 



dW/d<f> 



dej) ~ 50 



(6) 



gives us 4> ohs . 

Fifty or so e-foldings is the minimum amount of in- 
flation capable of giving our observed Universe, but cru- 
cial for the considerations of this paper is that typically 
one expects enormously more inflation. If one only im- 
poses the constraint that the initial energy density be 
below the Planck energy, most models allow very large 
amounts of inflation. However the early stages of such 
inflation may be dominated by quantum dynamics (the 
so-called stochastic inflation regime [12]). The maximum 
value of the inflaton field we consider, </> max , is given by 
the quantum-to-classical evolution transition 



H 
2^ 



H 



< 1 



(7) 



where the A represent the typical quantum and classical 
evolution per Hubble time. 

For simplicity we focus on a inflaton model with a 
power-law potential 



w(4>) = w 
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The former constraints lead to 
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III. 



QUINTESSENCE FIELD EVOLUTION 
DURING INFLATION 



The models are now fully defined and we can study 
the dynamics of each field, beginning with the classical 
evolution. 



A. Classical evolution 

We assume that the quintessence field is completely 
uncoupled from everything else, including the inflaton. 
The equations of motion are the Euler-Lagrange equa- 
tions for the inflaton and quintessence fields 



• dTU •• • dU 



and the Einstein equation 



H 2 



8tt 
3mp[ 



-U 2 + W(<f>) + l -Q 2 + V(Q) 



(13) 



(14) 
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Because of the very small value of Vo, the evolution of the 
Universe is determined by the inflaton. The quintessence 
field may dominate if its value is very small, but, as we 
will see this situation would end very quickly. Therefore, 
using the slow-roll approximation we have the solution 



0(N) 



9. 

mi 



aN 

47T 



(15) 



where N — log(a/a; n i). The subscript "ini" always means 
the value at some time Knowing the dynamics of the 
Universe we can now study the quintessence field. 

The classical behaviour of the quintessence field is 
rather straightforward. If the field is in a region where 
|dU/dQ| is small, its classical evolution will be highly 
suppressed by the friction arising from the inflationary 
expansion and to a good approximation the field will re- 
tain its initial value. Conversely, if |dU/dQ| is large, the 
field will quickly roll down until the potential becomes 
too flat. More precisely, using the slow-roll approxima- 
tion for the quintessence field we have 



Q 



PV 



2,H{t)ml x 



(16) 



Assuming H (t) constant (roughly correct for a few Hub- 
ble times) the solution is 



We define 



(t - iini) + Q 



0+2 
ini 



1/(0+2) 



Qr, 



P(P + 2) Vo 



'pi 



1/08+2) 



(17) 



(18) 



and see that the evolution during a Hubble time has the 
following approximate description: 



+ H- 1 ) 



— Qini 

— Qmin 



if Qini > Q n 
if Qini < Q 



(19) 



Therefore, after a few Hubble times the classical evolu- 
tion of the quintessence field is to remain constant. 



B. Quantum fluctuations 

Although the large friction from the Hubble expan- 
sion renders the classical evolution of the quintessence 
field negligible, the same is not necessarily true of the 
quantum fluctuations. Indeed, as we shall see, the effect 
of these in the quintessence potential is to drive typi- 
cal regions of the Universe to quite large values of the 
quintessence field, corresponding to low energy densities. 
This possibility was first mentioned in Ref. [13]. 

Let us now recall briefly the way we treat these fluc- 
tuations. Following Ref. [11], we choose the spatially- fiat 



slicing and we split the quintessence field into an unper- 
turbed part and a perturbation 



Q(x,t) = Q(t) + 5Q(x,t). 



(20) 



We quantize the perturbation and expand its Fourier 
components 



SQk{t) = w k {t)au + wl(t)a'_ 



k ' 



(21) 



where Sk is the annihilation operator. In the linear ap- 
proximation we have to solve 



w k + 3Hw k + [(k/a) 2 + m 2 ff ] w k 



0. 



(22) 



where m 2 ff = d 2 U/dQ 2 . If the field is effectively massless 
and H is approximately constant, the power spectrum of 
the fluctuations is given by 



V k (t) 



L 3 fc 3 
2tt 2 



\w k (t)\ 2 ~ 



k 
aH 



This means that before horizon-crossing V k (t) oc 1/ a 2 (t) 
and afterwards it becomes constant. Using 



(<5Q 2 (x,i)) 



v k (t)f, 



(24) 



and the fact that the fluctuations become classical after 
they cross the horizon, we find that the quintessence field 
receives quantum kicks whose size is estimated as H/2n 
per Hubble time. 

The effect of these quantum fluctuations can be stud- 
ied using the Fokker-Planck formalism [14], which allows 
one to follow the probability distribution f(Q,t) of the 
quintessence field during inflation. A simple derivation 
is as follows. We take discrete steps At = H~ 1 e during 
which there are a random jump of magnitude AQ qu = 
v /At/H~ 1 H/2Tr = \ r KeH/2ir due to quantum fluctua- 
tions and a classical step AQ c i = QAt = —sV'/3H 2 , 
where the prime denotes a derivative with respect to 
Q. During this time an interval dQ shrinks by a factor 
1 + sV"/3H 2 . We find the following equation: 



f(Q,t + At)dQ = -[f(Q-AQ ci 



AQ qul t) 



+ f(Q-AQ cl -AQ qu ,t) 
x (1 + eV"/3H 2 ) dQ . 



(25) 



Subtracting f(Q, t)dQ, dividing by e and taking the limit 
e — > leads to 



8f(Q,t) _H 3 (t) d 2 f(Q,t) 



dt 



8tt 2 



+ 



dQ 2 

d fdV(Q) 



3H(t) dQ V dQ 



f(Q,t) 



(26) 



On the right-hand side, the first term produces diffusion 
and the second one produces a drift. Note that unlike 



4 



the usual stochastic inflation situation, here H(t) is de- 
termined externally by the inflaton field evolution, and 
the fluctuations in Q do not back-react on the expansion 
rate. 

To be complete we need the boundary condition 



a=2 



lim 

Q—o 



H 3 (t)df(Q,t) , 0V(Q)/dQ 



8tt 2 OQ 



+ 



3H(t) 



-f(Q,t) 



-0, (27) 



the requirement that f(Q, t) > 0, and normalization 
Jo° f(Q, t)dQ = 1 which implies limQ^oo f(Q, t) = 0. 

From now on we focus on the inverse power-law po- 
tential for the quintessence. First of all, one may wonder 
whether the massless condition we used for the typical 
quantum jump is justified. Indeed, for values smaller 
than Qiim, where 



Q 



lim 



V03+2) 



(28) 



the quintessence field has an effective mass m e s which is 
no longer negligible compared to the Hubble rate. 1 For 
all wavelengths a/k bigger than l/m e g, Eq. (22) implies 



w k (t) oc a 



(29) 



and therefore the power spectrum of Pk(t) decreases as 
l/a 3 (t) leading to a smaller value of the quantum jump 
after horizon-crossing. Nevertheless, because the typical 
quantum jump H/2ir (typically 10 -5 mpi) is dramatically 
bigger than Qi; m (with 0=1, typically 10 _39 mpi), the 
massless condition is broken only during very short pe- 
riods of time (much less than a Hubble time) and so the 
massless approximation is an excellent one. 
The Fokker-Planck equation becomes 



df(Q,t) H 3 (t) d 2 f(Q,t) 



dt 




- 1 f(Q,t) 



(30) 



For Q > Q m i n we can drop the terms coming from the 
classical evolution since we have seen it is negligible. For 
Q < Qmin the effect of the potential becomes important 
and it acts as a wall preventing the quantum fluctuations 
driving the quintessence field to the origin. Therefore, to 
a good approximation we have to solve 



df(Q,t) _ H 3 (t) 2 f(Q,t) 



dt 



8tt 2 8Q 2 



(31) 



with the boundary condition df(Q w 0,t)/dQ = (no 
flux through the origin) and maintaining f(Q,t) > and 
Jo f{Qi t)dQ = 1. This is equivalent to a random walk 



1 For power-law potentials, this field value is almost the same as 
Qmim the value below which the classical evolution is important. 
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FIG. 2: Probability distributions for the mean quintessence 
value at the end of inflation for different initial values 0i n i, 
in the particular case of a = 2. We have also plotted the 
probability distribution obtained by averaging over <j) lVL \. 



with a wall. Starting with a (half) Gaussian distribution 
centred on Q = and with variance a 2 (t) , and assuming 
H(t) is constant during At, we have 



a 2 (t ini + At) = <j 2 ni + ^At 



More precisely, the equation to solve is 

H*(t) 



dt 



Att 2 



(32) 



(33) 



which has solution 



a 2 W 



+ 



16 Wo 



3a(a 



2)m P1 



(Cf-r +2 ). (34) 



After enough time the evolution is roughly independent 
of the initial condition crini and we can set it to zero. 

We can now ask what one expects the mean value of 
the quintessence field to be in our region of the Uni- 
verse. Because we are presently interested only in the 
mean value, we should only consider perturbations on 
scales larger than our present horizon, which are gen- 
erated between the initial inflaton value (/>; n i and ^> bs- 
Because of the presence of the wall, the net effect of the 
quantum fluctuations is to diffuse the quintessence field 
to larger values. The extent of this diffusion depends on 
how much inflation occurs before the last fifty e-foldings; 
if more inflation occurs then there are more 'steps' in the 
diffusion and also the early steps are larger. We have also 
computed the distribution obtained by averaging over all 
possible initial inflaton values between o b s and max as- 
suming a flat probability. 

Actually we are more interested in the probabil- 
ity distribution f q of q = log 10 (Q/mpi). We have 
f q (q,t) = f(Q(q),t)dQ/dq which has a peak at g pca k = 



W(<\>)=W Wm n ) 



V(Q)=V (Q/m n f 
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FIG. 3: The position of the peak of the probability distribu- 
tions of Q, as in Fig. 2, plotted as a function of a. 



log 10 (cr(i)/mpi). In Fig. 2 we show some distributions 
fq(q,t bs) m the particular case of a = 2; we show the 
probability distributions for two possible initial condi- 
tions for the inflaton, and then the probability distribu- 
tion averaged over a uniform initial distribution for </>;„; 
between <fi b s and </> max - Fig. 3 shows the position of the 
peak in the distributions as a function of a; as we can see, 
the result is roughly independent of a, and as expected 
the more inflation there is the further the distribution 
diffuses to large Q. 

We wish to know the quintessence value after inflation 
to assess when the solution begins tracking. In the slow- 
roll approximation, the solution for the tracker is 



Q(a)~ 



/ /?(/? + 2) y 

V247T/9fo(l + Wf) 



1/(0+2) 



3(l+io f )/(/3+2) 



(35) 



where pto is the present value of the energy density and 
Wf the pressure-density ratio of the dominant fluid (ei- 
ther radiation or non-relativistic matter). Some exam- 
ples of trackers are plotted in Fig. 4. After inflation the 
quintessence field remains constant as long as its energy 
density is lower than the tracker's. We can therefore find 
the probability distribution f z for log 10 ( z tr + 1), where 
Ztr is the redshift at which the quintessence reaches the 
tracking behaviour, by comparing the value of Q at the 
end of inflation with the tracking solution. Using Fig. 3 
and Fig. 4 we can estimate at which redshift the tracker 
is reached for some values of a and j3. We notice that 
the smaller [3 is, the later the tracker is reached for a 
given initial value of Q. In Fig. 5 we show some exam- 
ples of distributions for the redshift at which tracking 
begins. The discontinuity in the distribution comes from 
the fact that at about radiation-matter equality there is 
a transition between two trackers. 

The striking feature of this figure is how low the ex- 
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FIG. 4: Some examples of trackers for different values of (3. 
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FIG. 5: Probability distributions for the redshift at which the 
quintessence field reaches the tracker, for a = 2 and /3 = 1. 
The distributions are shown for two different initial values 
0ini, and averaged over </>i n i. The discontinuities arise across 
matter-radiation equality. 



pected redshift of tracking is. While some previous pa- 
pers have advocated equipartition of the quintessence en- 
ergy density as an initial condition [5] , leading to prompt 
tracking, we find that tracking is postponed until the late 
stages, and in particular well after nucleosynthesis. In- 
deed, the bulk of the probability is not only after nucle- 
osynthesis but after decoupling too; however there is only 
a small probability that tracking has not begun by the 
present, which would not lead to acceptable quintessence. 

Finally we must say that our results are valid for small 
initial values of the quintessence field before inflation (in 
fact formally zero). If this were not the case, the proba- 
bility distribution for the mean value of the quintessence 
field at the end of inflation will broaden to larger values 



6 



and the tracker will be reached even later or not at all. 



IV. CONCLUSIONS 

We have analyzed the quantum dynamics of the 
quintessence field during a period of early Universe in- 
flation. Due to quantum fluctuations, even if the initial 
value of Q in a certain region of the Universe is small, it is 
rapidly diffused to large field values and hence low energy 
density. We have found that typically tracking behaviour 
begins only at quite a late stage of evolution, well after 
nucleosynthesis and quite likely after decoupling too. 

Although we have discussed specific models of both in- 
flation and quintessence, we expect our results to be quite 
general; as far as inflation is concerned we need only the 
assumption that there are significantly more than fifty e- 
foldings in total and a standard value of H during infla- 



tion, while for other models of quintessence where the po- 
tential diverges at the origin the result should also be the 
same as the precise form of the potential is dynamically 
irrelevant. In particular, while we have only considered 
models where the tracking density during nucleosynthesis 
is negligible, these considerations may reinstate models 
whose tracking density would be unacceptable during nu- 
cleosynthesis. 
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